Noise-induced macroscopic bifurcations in globally-coupled chaotic units 
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Large populations of globally-coupled identical maps subjected to independent additive noise 
are shown to undergo qualitative changes as the features of the stochastic process are varied. We 
show that for strong coupling, the collective dynamics can be described in terms of a few effective 
macroscopic degrees of freedom, whose deterministic equations of motion are systematically derived 
through an order parameter expansion. 

PACS numbers: 05.45-a, 87.10. +c 



43 

o 

CD 



CO 



I 

C 

O 

o 



in 

in 
o 

o 

c3 



-a 
c 

o 
o 



X 



Dynamical units coupled all-to-all constitute relevant 
models of a wide range of physical and biological sys- 
tems such as Josephson junction arrays ll| , electrochem- 
ical oscillators 0] , interacting neurons |j| and yeast cells 
in a continuous-flow, stirred tank reactor 0. Globally- 
coupled dynamical systems are also widely used as pro- 
totypical models for addressing the collective behavior 
of systems with many degrees of freedom 0. In the 
limit of strong coupling of identical units, full synchro- 
nization usually occurs, i.e. the collective dynamics is 
an exact copy of the individual one. However, in the 
generic case where some diversity and/or noise is present 
at the microscopic level, macroscopic measurements do 
not directly reflect individual dynamical features. This 
simple observation is of large experimental importance. 
In the yeast cell experiment, for instance, single-cell mea- 
surements are not available so far. The idea is thus to 
extract information about metabolic processes at the in- 
dividual level from the observed collective oscillations in 
cell suspensions 0. 

In this Letter, we consider the influence of micro- 
scopic additive noise on the fully- synchronous regime of 
globally-coupled systems. This allows us to isolate the 
effect of noise from the possibly already complex collec- 
tive dynamics that systems of globally-coupled identi- 
cal and noiseless units may themselves exhibit 0. Re- 
stricting, for simplicity, our discussion to identical scalar 
maps, we show that noise plays a non-trivial role, "un- 
folding" the noiseless collective dynamics while progres- 
sively suppressing the direct connection between average 
and local evolution. We build a systematic expansion 
that accounts quantitatively for the collective dynamics 
observed, revealing that its essential features can be de- 
scribed in terms of few effective degrees of freedom and 
that its hierarchical structure is captured in progressively 
greater detail as the truncation order is increased. 

Our system of noisy and globally coupled identical one- 
dimensional maps is defined as follows: 

Xj ^(l-K)f( Xj )+K(f(x))+^(t) j = l,...N, (1) 

where Xj £ R denotes the state of the j-th popula- 



tion element, / : R — > R is a smooth function defining 
the dynamics of the uncoupled element, and ( f(x) ) := 
jf J2j=i f( x j) 1S tne average over the population. Every 
map is subjected to a white, but not necessarily Gaus- 
sian, noise of zero average and variance a 2 . 

Without noise and for strong enough coupling, the tra- 
jectories of all the population elements asymptotically 
coincide (full synchronization). The mean field X = (x) 
then evolves on a one-dimensional manifold of the high- 
dimensional phase space of the full system. 

Although noise smears any individual trajectory, for 
large population sizes these fluctuations average out and 
the mean field X evolves deterministically up to finite- 
size effects. Somewhat counter-intuitively, the macro- 
scopic dynamics is qualitatively different from that of the 
uncoupled map and remains apparently low-dimensional 
for any noise intensity. For chaotic logistic maps, increas- 
ing noise drives the mean field X through a cascade of 
period-doubling bifurcations which takes place at smaller 
and smaller a values as the coupling strength K is low- 
ered (Fig. n (d-f))- The difference between single ele- 
ment noisy evolution and regular average behaviour can 
be easly seen comparing the sharply peaked probability 
distribution functions (pdf ) of X with the broad, smooth 
pdfs of individual elements (Fig. [T](a-c)). 

The effect of noise is not always the "simplification" of 
the collective dynamics; rather, noise unravels the under- 
lying complexity of the phase space. This is particularly 
striking in the case of excitable units, such as: 



f(x) = («i+ 7 i J ) e"" 1 . (2) 

(For other examples of excitable maps, see Q.) While the 
origin is the unique global attractor, map <|2J) possesses 
an unstable fixed point, excitation across which triggers 
a chaotic transient. Increasing the noise intensity over a 
threshold a*, the mean field displays intermittent bursts, 
due to the crisis of a chaotic attractor which disappears, 
for stronger noise, via a period-doubling cascade (Fig. 
EJl. The onset of collective oscillations is nothing but 
an explicit mechanism for "noise-induced" and "coher- 
ence" resonance 8J: the single-element power spectrum 
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FIG. 1: System of N = 2 20 globally coupled logistic maps in 
the chaotic regime {f(x) = 1 — ax 2 , a — 1.57) subjected to 
uniformly-distributed noise of variance a 2 , (a-c): K = 0.7: 
pdfs for the mean field (solid line) and one single element 
(dash-dotted line; dotted line at even times, dashed line at 
odd times) for different noise intensities: (a) a 2 = 0.0025, (b) 
a 2 = 0.0324, and (c) a 2 = 0.09. (d): Bifurcation diagram for 
X as a function of a 2 for K = 1. (e): First period-doubling 
bifurcations for the mean field (o) and for the second-order 
reduced system (solid line, see text) in the (cr, K) plane, (f): 
Same as (d) but for the reduced system Eq.© with K = 0.7 
(the diagram of the full system almost coincides). The dotted 
lines in (d) and (f) indicate the location of the bifurcations 
reported in (e). 



displays an enhanced peak in the periodic regimes occur- 
ring at intermediate a values (ll). 

Let us now investigate analytically the apparent low- 
dimensional nature of the collective motion by trying to 
decouple the macroscopic effect of noise from the dy- 
namical "skeleton" provided by the microscopic map. 
We write the position of each population element in 
terms of the mean field and of its displacement from it: 
xj = X + ej . The uncoupled element equation can now 
be expanded as a power series around the mean field: 



oo 

f( Xj ) = f(X) + J2 -^fix) 



"3 ' 



(3) 




FIG. 2: Population of N = 2 22 excitable maps (Eq. © with 
a — 0.4, j3 = 1, 7 = 8, and K = 0.9) subjected to Gaussian 
noise of variance a 2 , (a): Bifurcation diagram for X vs a. 
(b): Critical value a* of the noise intensity over which the 
mean field displays large amplitude oscillations for the popu- 
lation (K = K = 0.8, o; K = 0.6, *) and for the order 
parameter reduction (solid line). 



where T> q f indicates the g-th derivative of the function /. 
Substituting Eq. into Eq. QJ, the macroscopic first 
iterate is obtained by averaging: 



oo 1 

X -</(£)> = /(*)+ E -7^/P0<^ 



(4) 



The mean field is then coupled to a full set of new order 
parameters: 



n q := 



q E N. 



(5) 



It is worth stressing that, at this point, we have taken the 
infinite-size limit, discarding the term (£) which scales 
like 1/y/N and plays the role of a macroscopic noise act- 
ing on the mean-field. This accounts for the mean field 
fluctuations obeying the law of large numbers. 
The evolution of the order parameters Q q can be com- 
puted using: 



oo 



ej -> (1 - K) \ - VPf(X) (e* - n p ) + 



P =i 



P 



and the fact that, in the limit N — > oo, the positions 
and the noise are uncorrelated variables and, hence, 
( h{X, e) £ 9 ) = ( h(X, e) ) ( ^ ) . Eq. © then yields: 



a, 




(1 - Kfrn^ 



(6) 



9=1 



where m q = (£, q ) is the g-th moment of the noise distri- 
bution. The macroscopic dynamical system defined by 
Eqs. J3J and © is still infinite-dimensional, like the full 
system, but it can be truncated to finite-order when as- 
suming that 1 — K is sufficiently small. This leads to 
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a hierarchy of finite-dimensional reduced systems which 
capture, as we show below, the properties of the collective 
dynamics in increasing detail. 

A first approximation of Eq. © consists in keeping 
only terms of zeroth order in (1 — K), so that any Q q — 



(a) 



The scalar equation for the mean field: 

oo 

X i ► f(X) + ]T - V\f{X) m« 

9=1 q ' 



(7) 



describes exactly the bifurcation diagram for K = 1 Fig. 
n (d). Here, the macroscopic effect of noise is a one- 
dimensional unfolding of the uncoupled element dynam- 
ics, whose nonlinearities become progressively more rel- 
evant as the noise intensity increases. In particular, if 
the map is polynomial, only a finite number of moments 
influence the mean-field dynamics. This suggests a cri- 
terium for dividing the noise distributions into classes 
according to the collective behaviour they generate. 
For logistic maps, Eq. Q takes the simple form: 



X 



aX 1 



(8) 



This map is easily rescaled to a logistic map in which 
a 2 plays the role of the usual nonlinearity parameter. 
Although Eq. JHJ provides a good approximation to the 
regimes of very strong coupling, it is independent of K. 
Hence, it cannot account for the changes with K of the 
period-doubling bifurcation lines (Fig. ^e)). Truncat- 
ing Eqs. (@J and JHJ to the second order in (1 — K) 
(the first order term vanishes for zero-average noise) the 
coupling strength appears as a parameter for the macro- 
scopic equations of motion. In this case, the higher-order 
moments obey a recurrence relation leading, when / is 
polynomial, to a two-dimensional map. For our example 
of logistic maps, this reads: 



X 

n 2 



1 - aX 2 

2 . (1- 



- a fl 2 
K) 2 a 2 



(9) 



X 



[m 4 - 6cr 4 + (4x 2 -n 2 + 6<j 2 )n 2 ] 



Figure W^e) shows that even for intermediate coupling 
strengths this equation accounts quantitatively for the 
main features of the mean field dynamics. 
Looking at the bifurcation cascade of Fig. ^(f), one could 
think that Eqs. 10 might be further reduced to a scalar 
map. This is not possible, however, as the mean field first 
return map is folded (Fig. and hence the macroscopic 
system must have more than one variable. Remarkably, 
our approximation at order (1 — K) 2 (Eq. © for logistic 
maps) reproduces this folded structure. For both the full 
system and the second-order approximation Eq. (jSJ), the 
distance dx between the first return map of X and the 
parabola identified by the zeroth order truncation scales, 
as expected, like (1 — K) 2 or a 2 (Fig.|3fc,d)) and vanishes 
on K — 1, where the zeroth-order truncation Eq. JSJ is 
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FIG. 3: Population of Fig. Q(a): First return map of X for 
the full system (dark dots), its zeroth-order approximation 
(solid line), and its second-order approximation (light dots) 
(K — 0.4 and a 2 = 0.001). (b): Q.& versus Q.2 for the same 
parameters (dark dots: full system, zeroth-order: triangle, 
second-order: solid line, fourth-order: light dots), (c) and (d) 
distances d x and dn defined in text (circles for the full system, 
dashed lines for the second-order approximation, solid line for 
fourth-order system), ((c): a 2 = 0.03; (d): K = 0.6). 



exact. The "thickness" of the transverse folded structure 
shown in Fig. [3Ja) scales similarly indicating that, 
apart from finite-size effects that may prevent us from 
observing it, the dimension of the macroscopic attractor 
becomes larger than that of the uncoupled element as 
soon as K < 1 and a > 0. 

Going one step further, it is easily seen that there are 
still features of the mean field dynamics that cannot be 
captured by Eq. ©. The macroscopic variables il q for 
the population satisfy only approximately the relations 
used for reducing the dimensionality of the truncated sys- 
tem. For instance, displays a complex dependence on 
O2 (Fig. Elb)), while it is constant at zeroth order, and 
linear at second-order. In turn, this complex structure is 
retrieved if Eqs. (@J and I© are truncated to fourth order 
in (1 — K). In the (f^, ^4) plane, the distances da sep- 
arating the second-order truncation from the full system 
or the fourth-order truncation are comparable and they 
scale like the expected (I — K) A and er 4 (Figs. 0Hc,d)). 

Extrapolating from the above results, we conclude 
that the presence of noise unfolds the fully-synchronized 
regime into an infinitely complex but hierarchically or- 
ganized structure whose main features are captured by 
the lowest orders of our approximation scheme. Never- 
theless, in the parameter region considered here, we ex- 
pect the collective chaos to remain embedded in a finite 
dimensional space. This is perhaps best seen when calcu- 
lating the Lyapunov exponents of the finite-dimensional 
maps approximating the full system. The maps possess 
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only one positive exponent in the chaotic region, whereas 
increasing the order of the approximation adds "new" , 
more negative exponents to those existing at previous 
order. These exponent values agree well with those nu- 
merically computed by simulating the Perron-Frobenius 
dynamics of the population pdf , so that we are con- 
fident that they do represent the collective dynamics of 
the full system. 

In general, when the local map / is not polynomial, the 
order parameter expansion of Eq. (J2J) contains an infinite 
number of terms even for K = 1, and it is then neces- 
sary to introduce the further closure assumption a <C 1 
to obtain closed form approximations. We can illustrate 
this for the excitable map example mentioned above. At 
second order in both (1 — K) and a, the approximation 
accounts well for the establishment of collective oscilla- 
tions. The linear stability analysis of the origin for the 
reduced system provides an analytical condition for the 
onset of the intermittent behavior, which is in good agree- 
ment with the simulations of the population (Fig. 0(b)). 

To summarize, the collective dynamics of infinitely- 
large populations of globally-coupled identical maps sub- 
jected to independent additive noise is deterministic and 
undergoes qualitative changes as the features of the 
stochastic process are varied. For strong coupling, the 
collective dynamics can be described in terms of a few ef- 
fective macroscopic degrees of freedom, whose equations 
of motion have been systematically derived through an 
order parameter expansion valid for any smooth map and 
noise distribution. Such a reduced system allows us to 
understand how the nonlinearities of the map interact 
with the specific features of the microscopic noise. At 
the experimental level, our work indicates that the ob- 
servation of unusual pdfs of local variables such a those 
shown in Fig.n mav be traced back to the presence of an 
underlying chaotic local dynamics. 

Our approach reveals that in a large region of parame- 
ter space the collective motion, in spite of being endowed 
with a single positive Lyapunov exponent, takes place on 
a hierarchically-organized complex attractor. We believe 
that this complexity is but another facet of the "anoma- 
lous scaling" properties of the cloud of points represent- 
ing the population in the local phase space which were 
uncovered recently by Teramae and Kuramoto [9(, both 
phenomena being rooted in the structure of the global 
phase space around the fully-synchronized solution. 

The method proposed here, together with Ref. flflj ]. 
allows to compare the macroscopic bifurcations induced 
by different sources of microscopic disorder, such as noise 
and parameter mismatch (quenched noise). In this last 
case, an order parameter expansion with other closure 
assumptions provides a different low-dimensional unfold- 
ing of the local dynamics, accounting for the observed 
diversity of the bifurcation scenarios jllj . 



Future work should consider how all these aspects 
evolve as one lowers the coupling strength towards the 
weak-coupling regimes considered by Shibata et al. 0] 
where collective chaos is high-dimensional. Our work 
can be extended to the case of continuous-time systems, 
where we expect a similar effect of microscopic noise |l3| . 
Finally, it would be most interesting to make the connec- 
tion with lattices or networks of locally-coupled maps, 
and the so-called non-trivial collective behavior they ex- 
hibit. There, as suggested in [J^j, the equivalent of the 
noise considered here lies in the finite-amplitude fluctua- 
tions introduced by the neighbors and the finite correla- 
tion length in the system. 
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